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BACKWARD STOCHASTIC DIFFERENTIAL EQUATIONS WITH 
RANDOM STOPPING TIME AND SINGULAR FINAL CONDITION 

By A. Popier 

Universite de Provence 

In this paper we are concerned with one-dimensional backward 
stochastic differential equations (BSDE in short) of the following 
type: 

Y t = £- I Y r \Y r \ q dr- [ ZrdB r , t>0, 

JtAr JtAr 

where r is a stopping time, q is a positive constant and £ is a T T - 
measurable random variable such that P(£ = +00) > 0. We study the 
link between these BSDE and the Dirichlet problem on a domain 
J)Cl and with boundary condition g, with g = +00 on a set of 
positive Lebesgue measure. 

We also extend our results for more general BSDE. 

Introduction. Let P) be a probability space, B = (Bt)t>o a Brow- 

nian motion defined on this space, with values in R d . (J-t)t>o is the standard 
filtration of the Brownian motion. Also given are r a {J-jj-stopping time, 
£ a real, ^--measurable random variable, called the final condition, and 
/ : Q x M+ x R x R d -> K the generator. 

We wish to find a progressively measurable solution (Y, Z), with values in 
M x R d , of the BSDE 

(1) ** = £+ r f(r,Y r ,Z r )dr- r Z r dB r , t>0. 

JtAT JtAT 

Such equations, in the nonlinear case, have been introduced by Pardoux 
and Peng in 1990 in [19], when r is replaced by a constant time T > 0. They 
gave the first existence and uniqueness result. Since then, BSDE have been 
studied with great interest (see the references in [18]). In particular, Peng 
[20] describes how the solution Y of (1) for an unbounded random terminal 
time is related to a semilinear elliptic PDE. Viscosity solutions for such 
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equations will be constructed by stochastic methods (see Theorem 8 below). 
This generalization of the Feynman-Kac formula is a reason for studying 
random terminal times. 

Let us recall the definition of a solution of (1) which can be found in [4]. 

Definition 1. A solution of the BSDE (1) is a pair {(Y t ,Z t ),t > 0} of 
progressively measurable processes with values in R x R rf such that, P-a.s.: 

• on the set {t > r}, Y t = £ and Z t = 0, 

• t h-> l t < T f(t, Y t , Z t ) belongs to £^(0, oo), t h-> Z t belongs to Lf oc (0, oo), 

• and for all < t < T, 



TAt 



YtAr = Ytat + f(r,Y r ,Z r ) dr 



Mr 



TAt 



Mr 



Zf dJ3 r . 



A solution is said to be an L p -solution for some p > 1 if, moreover, for some 
AgR, 

e( sup e pXt \Y t \ p + [ T e pXt \Y t \ p dt + f e pXt \Y t \ p ~ 2 \\Z t \\ 2 dt) < +oo. 
\0<t<T Jo Jo J 

We assume that the generator / : Vt x R + x R x R d — > R is such that: 

(HO) f(-,y,z) is progressively measurable, for all y,z; 
(HI) 3K > 0, such that a.s. Vt, y, z, z' , 

\f(t,y,z)-f(t,y,z')\<K\\z-z'\\; 

(H2) 3fi G R, such that a.s. Vi, y, y', z, 

(y-y'){f(t,y,z) -f{t,y',z)) < n\y-y'\ 2 ; 

(H3) y i— > /(i, y, z) is continuous, Vt, z, a.s. 

(H4) for all r > and all n G N*, </v(i) = sup|„|< r y, 0) - /(i,0,0)| be- 
longs to L 1 ((0, n) x fi). 

Now for some p > 1 we suppose that there exists A > v p = \i + 2 {p-i) > sucn 
that 



(H5) 

and 

(H6) 



E 



< +oo 



E 



|£|P + / e pXt \f(t,e-»o%,e-^%)\ p dt 



< +oo, 



where £ = e VpT ^, £ t = ^(£1 ^i) and r/ is predictable and such that 

,P/2"1 

e = E(0+/ %dB t , E 



l%| 2 ^ 



< oo. 



Let us recall Theorem 5.2 of [4]. 
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Theorem 1. Under the conditions (H0)-(H6), there exists a unique 
solution (Y, Z) of the BSDE (1), which, moreover, satisfies, for A > u p such 
that (H5) and (H6) hold: 

e( sup e pXt \Y t \P+ [ T e pXr \Y r \ p - 2 (\Y r \ 2 + \\Z r \\ 2 )dr) 

\0<t<T Jo J 

( 2 ) 

< cE + f e pXr \f(t, 0, 0) r dr) , 

/or some constant c = c(p, A, if, //). 

Remark 1. The previous theorem is a generalization of the result of 
Darling and Pardoux (Theorem 3.4 in [6]) or of Pardoux (Theorem 4.1 in 
[18]). In [6] or [18] the result is given in the case p = 2. Here we have expressed 
the theorem for the dimension one (£ and Yf belong to M). But it is still true 
in higher dimensions (see [4]; the product in (H2) must be replaced by the 
scalar product in K m ). 

Note that if / is a Lipschitz function, the condition (H2) holds. 



From now and in the rest of the paper we are concerned with the BSDE 

(3) Y t = £- [ T Y r \Y r \ q dr - [ Z r dB r with q > 0. 

JthT JthT 

Here the function / is deterministic and equal to 

f{t,y,z) = -y\y\ q . 

f satisfies all conditions (H0)-(H4) of Theorem 1, with K = \i = (which 
implies u p = for all p > 1). Indeed, / is a nonincreasing function, thereby, 

-(y -y')(y\y\ q -y'\y'\ q ) <o. 

Since /(t, 0,0) = 0, (H5) is always satisfied. 

The stopping time r is defined as follows. Let D be an open bounded 
subset of whose boundary is at least of class C 2 (see [12] for the definition 
of a regular boundary). For all x £ M. d , let X x denote the solution of the SDE: 

(4) Xf = x+ f 'b{X*)dr+ f a(X?)dB r for t > 0. 

Jo Jo 

The functions b and a are defined on M. d , with values respectively in M. d and 
]R rfxd , and are measurable such that: 

• Lipschitz condition: there exists K > such that 

(L) V(x,y) £R d xR d \\a(x) - a(y)\\ < K\x - y\; 



4 



A. POPIER 



• Boundedness condition: 

(B) VxeR d \b(x)\ + \\a(x)\\ <K; 

• Uniform ellipticity: there exists a constant a > such that 
(E) VxeR d aa*{x)>ald. 

In the rest of this paper (L), (B) and (E) are supposed to be satisfied. 
Under these assumptions, from a result of Yu Veretennikov [24] and [25], 
equation (4) has a unique strong solution X x . For each x 6 D, we define the 
stopping time 

(5) t = t x = inf {t > 0, Xf i D}. 

Our stopping time satisfies the following two properties. Since D is bounded 
and since the conditions (B) and (E) hold 

(CI) every point x € dD is regular. 

In particular, if a; € dD, t x = a.s. (see [3], Corollary 3.2). This assump- 
tion (CI) is important to define a singular solution (see Definition 2 below). 
Moreover, since (L), (B) and (E) hold, we have the following result (see [21], 
Theorem 2.1 and [18], Remark 5.6): for all x S D, t x < +oo a.s. and there 
exists (3 > such that 

(C2) supE^*) <oo. 

x£D 

This property will be used to construct solutions of the BSDE (3) for 
bounded terminal conditions £ (see Proposition 2). 

From the papers [6, 18] or [20], we know that the BSDE (3) with terminal 
time equal to r = t x and final data equal to £ = h(X^ ) is associated with 
the following elliptic PDE with Dirichlet condition h: 

— Cu + u\u\ q = on D, 

(6) 

u = h on dD; 

where C is the second order partial differential operator: for all ip G CQ(M. d ), 

(7) Vx G M. d £tp(x) = ±Trace(aa*(x)D 2 (p(x)) + b(x)V<p(x). 

In the rest of this paper V and D 2 will denote respectively the gradient and 
the Hessian matrix. If (Y X ,Z X ) denotes the solution of the BSDE (3) with 
terminal data h(X® ), the connection is given by the formula 

u(x)=Y x . 

Le Gall [13] succeeded in describing all solutions of the equation Au = u 2 
in the unit disk D in R 2 by a purely probabilistic method. He established 
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a 1-1 correspondence between all solutions and all pairs (T,u), where T 
is a closed subset of dD and v is a Radon measure on dD \ T. The set 
r is the set of singular points of dD where the solution explodes badly: 
roughly speaking, near points of T, the solution behaves like the inverse of 
the squared distance to the boundary. The measure v can be interpreted as 
the "boundary value" of u on dD\T. The solution corresponding to (T, v) is 
expressed in terms of the Brownian snake (a path- valued Markov process). 
In [14] the results announced in [13] are proved in detail and are extended 
to a general smooth domain in M 2 . 

The pair (T, u) is called the boundary trace for positive solution of the 
PDE (6). The definition of boundary trace in general was provided by Mar- 
cus and Veron [15] who showed by analytic methods that every positive 
solution of (6) possesses a unique trace. The trace can be described by a 
(possibly unbounded) positive regular Borel measure v on dD. The corre- 
spondence between (T, v) and u is given by 



for every Borel subset A of dD. 

The corresponding boundary value problem is presented in [15] in the 
subcritical case < q < 2/(d — 1) and in [16] in the supercritical case q > 
2/(d — 1). In the subcritical case, for every pair (T,v), the problem has a 
unique solution. Remark that in [13] and [14], q = 1 and d = 2, that is, the 
subcritical case is studied: q = 1 < 2/(2 — 1) = 2/{d — 1). In the supercritical 
case Marcus and Veron derive necessary and sufficient conditions for the 
existence of a maximal solution. Similar conditions were obtained by Dynkin 
and Kuznetsov [8] for q < 1. Their method relies on probabilistic techniques 
and is not extendable to q > 1, because the main tool is the g-superdiffusion 
which is not defined for q > 1 . 

The object of the present paper is to give a probabilistic representation of 
the solution of the PDE (6) in terms of the solution of the related BSDE (3). 
In general, a solution of the PDE has a "blow-up" set T. Therefore, the final 
data £ of the BSDE must be allowed to be infinite with positive probability 
and the set {£ = +00} corresponds to T. Hence, our first problem is to find 
a solution of (3) when £ is infinite with positive probability, which implies, 
in particular, that (H6) is not satisfied. 

Note that there are some differences between our work and the results of 
Le Gall or Dynkin and Kuznetsov. With the superprocesses (see [14] or [8]), 
it should be assumed that q < 1. In our case there is no restriction on q > 0. 

Moreover, the Dirichlet boundary condition for the PDE (6) is not taken in 
the same sense in the two approaches. With the notion of the boundary trace 
(see [8, 14, 15] and [16]), there always exists a maximal positive solution; if 
q < 2/(d — 1), this solution is unique, and if q > 2/{d — 1), the problem (6) 




if AQ{dD\T) 

if 
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may possess more than one positive solution. More precisely, assume that 
D is the unit ball in that q > 2/(d — 1) and denote by p^ the Borel 
measure on dD which assigns the value +00 to every nonempty set. Then 
for every e > 0, there exists a positive solution of (6) such that u(0) < e and 
the trace of u is p^ (see Proposition 5.1 of [16]). 

In our case the Dirichlet condition in (6) is taken in the viscosity sense 
(see Definition 4 in Section 5). The results are rather different: there exists 
a minimal positive viscosity solution. But we are unable to give conditions 
to ensure uniqueness of the solution. 

Main results. In the first section we will prove an a priori estimate which 
is a probabilistic generalization of the Keller-Osserman inequality. 
In Section 2 and in the rest of the paper we assume 

f > a.s. 

and we allow £ to be infinite with positive probability: P(£ = +00) > 0. We 
must modify Definition 1 of a BSDE when £ does not satisfy the condition 
(H6). 

In the rest of this paper p denotes the distance from the boundary of D. 
For x £ D, for all positive 77, let us define the stopping time 

(8) r*=mf{t>0,p(X?)< V }. 



Remark 2. For x 6 D, < t x a.s. and if x 6 D, when r] goes to 0, 
converges to t x a.s. When x € dD, for all 77 > 0, = t x = a.s., because 
every point x € dD is regular [condition (CI)]. 

Therefore, we suppose x to be in D and for convenience, we omit the 
variable x. 

Definition 2. For an ^--measurable £ such that P(£ > 0) = 1 and 
P(£ = 00) > 0, the process (Y, Z) is a solution of the BSDE (3) if: 

(Dl) for all 77 > and all T > 0, 

/ /-TAt,, \ 

E sup \Y tA T | 2 + / \Z r \ 2 dB r <+oo; 

\0<t<T J0 J 

(D2) P-a.s. for all < t < T and all 77 > 0, 

*W„ = Y TATv - Y r \Y r \«dr- / " Z r dB r ■ 
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(D3) on the set {t > r}, Yt = £ and = 0, and P-a.s., 

lim Y tA T = €- 

t—t+OO 

We first construct a process {(It, ^t);i > 0} satisfying the conditions (Dl) 
and (D2) of the previous definition. (Y, Z) is the limit of the sequence of 
processes (Y n ,Z n ), solution (in the sense of Definition 1) of the BSDE (3) 
with terminal condition £ A n. From the first section we already know that 
there exists a constant C such that 

Vt>0 p 2/q (X tAr )Y t <C. 

Moreover, we prove the following: 

Proposition 1. On {£ = +00} the explosion rate ofY is in the order 
of p~ 2 l q (Xt/\ T ): there exists a positive constant C depending on D, q, the 
bound on b and a in (B) and on the constant a in (E), such that 

liminf p 2 /i(X tAr )Y tAT > C a.s. on {£ = +00}. 

Without other assumption on £, we cannot prove that (Y, Z) satisfies the 
condition (D3) of Definition 2. 

In Section 3 we first prove that P-a.s. the limit of Yt Ar as t goes to +00 
exists and 

lim Y tAT > £. 

t— >+oo 

Then we add some assumptions on £ and on the diffusion X to insure that 
the condition (D3) holds. We prove the following: 

Theorem 2. Under the assumptions: 

• the terminal data £ satisfies 

(Ai) Z = g{Xr), 

where g : M. d — > M + is a function such that = {g = +00} P\dD is a closed 
set; 

• on M. d \ Fqo, g is locally bounded, that is, for all compact set K, C M. d \ F^, 
(A2) gl K eL™{M. d ). 

• the boundary dD belongs to C 3 ; 

the process Y is continuous, that is, linif^ +00 Yt Ar = £ P-a.s. 
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In the next section we prove if there exists a solution (Y, Z) of the BSDE 
(3) in the sense of Definition 2, then Y > Y. Therefore, if the process (Y, Z) 
is a solution (e.g., if the assumptions of Theorem 2 hold), it is the minimal 
solution. 

In the last section we show the connection between the BSDE (3) with 
terminal condition g(X®) and the PDE (6) with Dirichlet condition g. The 
assumptions of Theorem 2 hold. In the previous sections we have defined a 
process {(Y t x , Zf);t >} which is a solution of the BSDE (3) with terminal 
data g(X% x ). Next we define 

u(x)=Y*. 

The main result follows: 

Theorem 3. Under the assumptions of Theorem 2, u is a viscosity 
solution of the PDE (6) with Dirichlet condition g. 

Here we do not suppose that a viscosity solution is continuous. But under 
some stronger assumptions on the operator C, we also give some regularity 
properties of the solution it. We also prove that u is the minimal solution. 

In the last section we will see that these results are still true with more 
general generators /. 

Theorem 4. Assume that f is a nonincreasing and C 1 function with 
f(0) = 0, and such that there exists q > 0, k > s.t., 

(9) Vy>0 f(y)<- K y 1+q . 

If £ is a nonnegative random variable, with P(£ = +oo) > 0, and such that 
the assumptions of Theorem 2 hold, then there exists a process (Y,Z), solu- 
tion of the BSDE 

(10) Y t = £+ f f(Y r ) dr - [ T Z r dB r 

JtAr JtAr 

[in the sense of Definition 2, with f instead of y i— ► —y\y\ q in (D2)]. 

Moreover the conclusion of Theorem 3 is still true: there exists a minimal 
viscosity solution for the PDE 



(11) 



Cu + f(u) = onD, 
u = g on dD. 
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Important remark on the condition (E). The condition (E) can be re- 
laxed. In the rest of the paper we can also work with the assumptions (L), 
(B) and we add the following condition: b is continuous and satisfies the 
monotonicity condition: there exists \i € M. such that 

(M) V(x, y) £ R d x R d (x - y\b(x) - b(y)) < fi\x - y\ 2 ; 

here (•[•) denotes the scalar product in M. d . Under these assumptions equation 
(4) has a unique strong solution X x . For each x € D, we define the stopping 
time 

T = T x = M{t>0, X? £D). 

We also assume that the conditions (CI) and (C2) hold. 

Under the assumptions (M), (L), (B), (CI) and (C2), the results which 
may be false are in Section 2, Proposition 1, and in Section 5, Propositions 
11 and 12. We are unable to control the explosion rate of Y (see Remark 5), 
nor to prove that the viscosity solution u is continuous on D without the 
ellipticity condition. 

In the rest of the paper all results (except maybe Propositions 1, 11 
and 12) could be proved without the condition (E). Indeed, we use this 
assumption only in the proofs of Propositions 4 and 8, in order to control 
the Green function G(x, •) associated to the process X x killed at t x . Under 
(E), this function G(x, •) is continuous on D except at the point x, and is 
integrable on D. This assumption on G can replace (E) (see, e.g., [21] for 
more details on G). 

1. An a priori estimate. Let (Y, Z) be the solution of the BSDE (3) with 
terminal data ^ such that the hypothesis (H6) holds. We will need an a priori 
inequality in order to control Yjat for t € [0, +oo[. The idea comes from the 
Keller-Osserman inequality which is true for any open set D ([11] and [17]). 
Denote by p the distance to the boundary of D C M. d . 

Theorem 5 (Keller-Osserman). There exists a positive constant C = 
C(q,d) such that if u is any C 2 (D) solution of 

-Au + u\u\ g = inD, 

then for all x € D, 

KX)! -^- 

Recall that in our case D is supposed to be bounded and dD £ C 2 . The 
process X x is the solution of (4), and the stopping time r x is defined by (5). 
We will prove the following: 
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Theorem 6 (A priori estimate). There exists a constant C [depending 
on the open set D, on q and on the bound in (B) of b and a] such that for 
every x £ D and every solution (Y, Z) of the BSDE (3) with terminal time 
t x and terminal data £ such that (H6) holds, we have 

C 

(12) Vi>0 \Y t \< 



We define the signed distance d 



, _ J dist(x,dD) = p(x), if xeD, 
[X) ~ \-dist(x,dD), if xeR d \D. 



For jj, > 0, let 



r^{xet d ;|rf(x)|</i}. 



The following lemma (see [9], Lemma 14.16) relates the smoothness of the 
distance function d in to that of the boundary dD. 

Lemma 1. Let D be bounded and 3D £ C k for k>2. Then there exists 
a positive constant fi depending on D such that d£ C fc (r M ). 

Proof of Theorem 6. Recall that D is an open bounded subset of M. d 
with dD £ C 2 . From the previous lemma we already know that there exists 
a positive constant fi such that on T^, the signed distance function d belongs 
to C 2 . And d = p is continuous on D. There exists a positive constant R 
(depending only on D) such that for all x £ D, < d{x) = p(x) < R. Let 
$ £ C°°(R d ; [0, 1]) such that $ is equal to 1 on l d \ T M and is equal to on 

r /V2- 

For < e < 1 and C > 0, we define a function ^ £ £ C 2 (R d ; R+) such that 
on D, 

* ° 



■ E 



[{l-$)p + R$ + E\ 2 /i' 



Such a function exists because (1 — <£)p + RQ + e > e on D. Remark that if 
x ££>, 

m £ (x) < 



g 



p(x) 2 / 

We denote by 6 £ the function (1 - <£)p + R® + e, that is, on Z), ^ 



2 ^ 9 . We apply the ltd formula to ^ £ {Xf hTx ), where x £ D. For conve- 
nience, we fix e > and x £ D and we omit the index e and x. For all 
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0<t<T, 



pTAt 

*(X tAr ) = V(X TAT ) - / ^> 1+q {X r ) dr 

JtAr 



(13) 



TAt 

tAr 
TAt 

Mr 



V^(X r )a(X r )dB r 



[V^(X r )b{X r ) 



+ \ Trace (aa*(X r )D 2 ^(X r )) - ^(X r ) 1+q ] dr. 



Now 



q,l+q _ £-1+90-2/9-2 
9 * 



dxi dxj q 



dxidxj q \q 
Therefore, 



dxi dxj 



(V*)6 + - Trace(ao-*L> 2 ^) - * 1+< ? 



(14) 



_ cg -y q -2 



9f) 

c q + —(ve)b 



1 /2 



- - + 1 llffVeir + - Trace (aa*D 2 9) 
q\q J q 

and b, a, 9, V# and D 2 9 are bounded on D. So we can choose the constant 
C such that for all x 6 D 



(15) 



+ lj ||ff(x)Ve(x) 



+ ^Mli:ace(<7(x)a*(x)£> 2 fl(x)) > 0. 
9 



The constant C depends only on D, on g and on the bound in (B) of b and 
a. We have obtained for all < t < T, 

rTAr 

#(X tAr ) = V(X TAT ) - / W(X r )a(X r ) dB r 

JtAr 

(■TAt (TAt 

- / ^{X r ) 1+q dr+ / U r dr; 

JtAr JtAr 
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with U a nonnegative adapted process, and on {t > r}, 

C 



£ 2/q- 

If (y, Z) is the solution of the BSDE (3) with a final condition £ in L°°(ft, .F r) p) 
(see Remark 4) , we can find < e < 1 such that 

Moreover, the Tanaka formula (see [10]) leads to, for all < t < T, 

i-TAt i-TAt 

\Y tAT \ = \Y TAT \ - / sign(Y r )Y r \Y r \ q dr - / sign(Y r )Z r dB r 

JtAr JtAr 

+ 2(A tAr (0)-A TAr (0)) 

pTAr i-TAt 

= \Ytat\- \Y r \ l+q dr- sign(Y r )Z r dB r 

JtAr JtAr 

+ 2(A t (0)-Ar(0)), 

where A is a local time of Y. Thus, At At (0) > At AT (0) a.s. 

By a comparison theorem (Corollary 4.4.2 in [6]), we have a.s. 

Vt>0 \Y t \ < y e (X tAT ) < - ° —. . 

tATj ~ p{X t AT?l q 

By a density argument, it is clear that, if (Y, Z) is the solution of the BSDE 
(3) with a final condition £ satisfying (H6), then 

C 



Vt>0 \Y t \< 



p(x t Ar?i q ' a 



Remark 3. The constant C depends only on the bound in (B), on D 
and q. Moreover, in the special case where D is a ball and where the drift in 
the SDE is equal to 0, the constant C depends only on q and on the bound 
of a and not on the center nor the radius of the ball. For example, if X is 
the Brownian motion, (12) is true for any C such that 

CO*(> + l) + ± 

Proof. In the case of a ball, we can give a slightly different proof be- 
cause we have an explicit expression for the distance function. We will as- 
sume that D is the ball centered at y and with radius R. In this case the 
function p is equal to 

. „ R 2 — \x — y\ 2 , . 
p(x) =R-\x-y\< = 6{x), 
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if x is in the ball. The function 9 is not of class C 2 on the whole space. We 
modify 9 in order to have a C 2 function. For all < e < R 2 , on the ball G, 
9 e will be equal to 

a r \ R 2 + £-\x-y? 



R 

and on the whole space K rf , 9 £ is positive and of class C 2 . Remark that on 
D, 9 e is greater than 9. Now we consider the function 

* e (x) = 



and like in the proof of Theorem 6, we prove that there exists a constant C 
such that the inequality (15) with 6 = holds. But now for i£D, 

Ve E {x) = ~{x-y) \\a{x)V9 E {x)\\ 2 = ^\\a{x){x-y)\\ 2 

<^\x-y?<AK 2 - 
D 2 9 £ (x) = -^ld =► Trace(cj(x)cr*(a;)L> 2 6 , e (a;)) 



i? 



It suffices to choose C such that 



2 2K 
-Trace(a(a;)cj*(x)) > - — . 



Q\Q J q 

in order to have (15). If X is the Brownian motion U, that is, tr = Id, C q > 
|( - + 1) + Here C depends only on the dimension d and on g, like in the 
Kell er-Osserman theorem. □ 



2. Approximation. We first prove a technical result which gives a suf- 
ficient condition on £ to insure existence and uniqueness of the solution of 
the BSDE (3). In our case for all p > 1, v v = and f(t, 0, 0) = 0. 



Proposition 2. Under the condition (C2) on £/ie ./irsi emf izme r o/ 
i/ie diffusion X , let £ be an J- T -measurable random variable such that £ £ L' r 
lozi/i r > 2(1 + q). Hence, there exists p > 1 and A > s-uc/t that the condition 
(H6) is satisfied: 



(H6) 



E 



» pXt \E(Z\F t )\ p{1+q) dt 



< +oo. 
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Proof. With a > 1, 7 > 1 such that l/a + l/j = 1, the Holder inequal- 
ity leads to 

E(e pXr \£\ p ) < [E(e apAr )] 1/Q [E(|£| 7P )] 1/7 . 

If £ £ L r with r > 1, there exists 7 > 1 and p > 1 such that -fp <r. From 
(C2), we can choose A > such that Xap < /3. 

For the rest of the condition (H6), we have f(y) = — y\y\ p and thus, 



E 



e pXt \f(E^^\ p dt 



<E 



<E 



e pA 'E^(|Cr (1+9) )dt 

p\r _ 1 



Xp 



sup E^(|e| p(1+9) ) 

te[0,r] 



1 



e sup E-^d^ 1 ^)) 

te[o,r] 



1/7 



From the Burkholder-Davis-Gundy inequality we obtain 

E sup E^(|£p p ( 1+g )) < C[E|e| 27p{1+9) ] 1/2 . 

te[o,r] 

If r > 2(1 + q), we can choose 7 > 1, p > 1 and A > sufficiently small such 
that 27p(l + q) < r and aXp < f3. □ 

Remark 4. From the previous proposition, if £ € U for some r > 2(1 + 
g), there exists p > 1 and A > such that the condition (H6) is satisfied. 
From Theorem 1, there exists a unique L p -solution (Y,Z) of the BSDE (3) 
which satisfies the estimate (2). 

For terminal data £, if (Y, Z) is a L p -solution for some p > 1, then (Y, Z) 
is also a L p -solution for all 1 < p' < p. Therefore, if (Y 7 , Z') is a L p -solution 
for some p' <p, then we can easily prove that (Y, Z) = (Y ; , Z'). 

If £ € L°° , from the proof of Proposition 2, (H6) holds for every p > 1 and 
A = /?/p- 



We are interested in the case where £ is a nonnegative random variable 
with this new assumption: 

P(£ = +oo)>0. 

We still assume that the conditions (L), (B) and (E) hold, that r = t x is the 
exit time of the diffusion X x from the set D, and that r satisfies (CI) and 
(C2). From Remark 2, we can suppose x to be in D and for convenience, we 
omit the variable x. We suppose that £ is ^--measurable. 

Now for each n £ N* , let £ n = £ A n be our final condition. With Remark 
4 we obtain the following: 
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Lemma 2. There exists a unique solution (Y n ,Z n ) (in the sense of Def- 
inition 1) of the BSDE (3) with terminal data £An. 

From a comparison theorem (see Corollary 4.4.2 in [6]), we have a.s. 
Vt>0, n<m 0<Y t n <Y t m . 
Define the progressively measurable process Y by 
(16) Y t = lim y t n V< > 0. 



Proposition 3. The sequence (Z n ) ng N* converges also to a process Z 
and (Y,Z) satisfies the assumptions (Dl) and (D2) of Definition 2. 

Proof. From (16), we already know that (y n ) n£ N* converges to Y. 
From the It 6 formula and the Burkholder-Davis-Gundy inequality, there 
exists a constant K such that for all rj > 0, n>m and < s, 



El sup \Y- -Y™ f)+E l SATV \\Z?-Z™fdr 



te[o,s] 



■ t ATjj ± t/\T v I 



<KM(\Y"-Y™ T f) 



See (8) for the definition of r^. But with the inequality (12) 



Y n ° ° 



r " ~ p(.X sATv )Vi ~ rj 2 /*' 

and with the Lebesgue theorem, we conclude that (Y™ T , T ) n converges 

to (Y ATv ,Z. ATv ) in L 2 (Q;C(R + ;R+)) x L 2 (0 x R+) and (Y% Tti ) converges 
uniformly to Y At . 

Hence, (Y,Z) satisfies the following equation: for all < t < s, for all 
77 >0, 



rsAT v rsAT v 

Y tATv =Y sATv - (Y r ) 1+9 dr- Z r dB r . 

J f At,, J tAr v 

From this equation, with the Ito formula and the estimate (12), we deduce 
that 

e( sup \Y tAT \ 2 ) +E [ SATV \\Z r \\ 2 dr < K~E\Y sAt | 2 
\te\o,s] J Jo 



K K 

< —T-, < 



Therefore, (Y, Z) satisfies the conditions (Dl) and (D2) of Definition 2. □ 
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From Definition 1, we also have that on the set {t > r}, Yj = £ and Zt = 0. 
With the monotonicity of the sequence Y n , we can conclude that 

liminf Y-t^ T > £ a.s. 

It remains to show the converse inequality, 

limsupY tAr < £, 

t—*+oo 

to have the last condition (D3). Without more assumptions on £, we cannot 
prove (D3). But we are able to give some other estimates on Y and Z. 

Proposition 4. For all e > 1, there exists K such that 
E 



1 [ \\Z r \\ 2 p{X r fl q+e dr<K. 
Jo 

Proof. We use again the notations T^, the function 9 = 9q = (1 — &)p + 
as in the proof of Theorem 6, and for r] < p. Recall that 9 € C 2 (D; M) 
and on D, 9 = (1 - $)p + i?$ > d > 0. Of course, x h-> |6»(x)| 4/9+e is not 
in C 2 (]R rf ), but this function belongs to C 2 (D \ F„) and we can define this 
function on the rest of (R d \ D) U T v in order to have the required regularity. 
The ltd formula leads to 
(Y t n ATv f9(X tAT y<> + c 

= {Y Q n ) 2 9(X ) 4 / q+£ + / \\Z?\\ 2 e(X r )V q+£ dr 
Jo 

/•(At, rtAT v 

+ 2 / (Y r n ) 2+q 9(X r )^ q+£ dr + 2 Y r n 9(X r )^ q+£ Z? dB r 
Jo Jo 

/A \ rtAr v 

+ \q +£ ) Jo ( y ") 2 ^(^r) 4/9+£ - 1 V0(X r )(6(X r )dr + a(X r )dB r ) 
rtAr n r/4 \ 

■ J Q (Y r n ) 2 [ [~ + e - 1 J 9(X r ) i / q+£ - 2 \\a(X r )V9(X r )\\ 2 
+ 9{X r fl q+£ - 1 Ttace(aa*(X r )D 2 f(X r )) 



(4/q + e) rtAT " 



f A \ rt/\T n 



Then 



ft At. 

E 



/ ||Z r n || 2 0(X r ) 4/9+e (ir 
Jo 

/ 4 \ At,, 

+ 2^- + ejEj Y r n 9(X r ) A / q+£ ~ l Z?V9(X r )o-(X r ) dr 
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is bounded from above by 

E((Y t n ATv fe(x tATv )^ 

'4 



+ ejE {Y?f9{X r ) 4 l q+£ - l V6{X r )b{X r )dr 

(17) _i(_ + e )E/ (Y r n ) 2 e(X r ) A / q+£ ^ 1 Tiace(aa*(X r )D 2 9(X r ))dr 
2 V q J Jo 

+ e) Q +e - ljE f^ Tr ' \Y r n ) 2 f(X r )^9(X r y- 2 



x \\a(X r )V0(X r )f dr. 

In the proof of Theorem 6, we have obtained that there exists some con- 
stant C such that for all n G N* and for all t > 0, 

(18) (l^ r ) 2 0(AW) 4/9 < C. 

Moreover 6 and cr are bounded [assumption (B)], and V6 and D 2 9 are also 
bounded on D. Thus, the right-hand side of (17) is bounded by 

K^I + eJ^ 8 £ ~ l (X r ) dr + E jT 9 £ ~ 2 (X r ) dr^j . 

We denote by p(t, x, y) the density of P x (X t 6 dy; r>t). P x means that the 
diffusion process X starts from x € D at time 0. Then 

E e £ -\X r )dr= / e £ - 1 (y)p(r,x,y)dydr 
Jo Jo Jd 

r~ 1 (y)G(x,y)dy, 

D 

where G is the Green function associated to the process X killed at time r 
(see [21], Section 4.2, Theorem 2.5). 

We claim that the last integral is finite. Indeed, if B{x,v) is the ball 
centered at x with radius v > 0, since x S D, we can find v > such that 
B(x,v) C D and B(x,v)nT u = 0. We denote by U the set D\(B(x,i/)UT u ). 
On U, f £ ~ 1 G(x,-) is a continuous function and is bounded by K. On T v 
[resp. on _B(x,z/)], G(x,-) (resp. 6> e ~ 1 ) is continuous and bounded by K. 
On the boundary of D, O 6 ^ 1 is singular if e < 1. Recall the definition of 
: = (1 — ip)p + Rip. Hence, is equivalent to p at the boundary and if e > 0, 
f £ ~ l is integrable on D. G has a singularity when y — > x, but with Theorem 
2.8 and Exercise 4.16 in [21], this singularity is integrable. Therefore, we 
split the integral into three terms: 



El 9 £ -\X r )dr= l D 9 £ - 1 (y)G(x,y)dy 
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3E-1/ 



B(x,v) 



+ 



L (y)G(x,y)dy + 



(y)G(x,y)dy 



QE-1 



(y)G(x,y)dy 



<K [ G(x,y)dy + K [ 6 £ ' x {y) dy + K Vol(D) 

JB(x,u) JYy 

< +OC. 

From the second integral, the same arguments show that if e > 1, 
(19) E [ T 6 £ ~ 2 {X r )dr <+oo. 

Therefore, the right-hand side of (17) is bounded by a constant K which 
does not depend on 77, n and t. And using the Cauchy-Schwarz inequality, 



E 



/ Y™0 (X r ) / q+e ~ Z^V0{X r )a (X r ) dr 
Jo 

< f E ^ ||Z r n || 2 fl(X r ) 4/<?+£ drJ 

E J* ATV (Y?) 2 8(X r )*/ q+£ - 2 \\V8(X r )a(X r )\\ 2 dr 



1/2 



But since V# and <r are bounded and since (18) and (19) hold, if e > 1, 



E 



(At„ 



(Y r n ) 2 8(X r ) 4 / q+£ - 2 \\V8(X r )a(X r )\\ 2 dr < K. 
Inequality (17) can be written as follows: A n + B n < C n with 

ftAT v 

0<A n = E \\Z?\\ 2 0(X r ) 4 /* +e dr, 

Jo 



\B n \ = 2 -+e 



E 



rtAT v 

/ y/wx^/^-^wfXrWXr) dr 

Jo 



< ka 1 J 2 



and |C„ I < K. Thus, for all n G N* and for all t > 0, 



E 



/ ||Z r n |j 2 &(X r ) 4 ^ +£ dr <K, 
JO 



which implies, by Fatou's lemma, 



E 



\Z r \\ 2 8(X r ) 4 /i +e dr < K. 



Since 8 > p on D, we obtain the announced result. □ 
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The condition e > 1 is required in order to insure that 

E f 9 £ - 2 {X r )dr < oo. 
Jo 

But this integral is equal to J D 6 e ~ 2 (y)G(x, y) dy, where G(x, y) is the Green 
function associated with the process X x killed at r. And if, for example, the 
infinitesimal generator of the diffusion X is self-adjoint in L 2 {M. d ), that is, 
£= (l/2)div(<7«7*V), then G(x,y) < Kp(y) (see [7], Theorem 9.5) and the 
previous integral is finite for any e > 0. Hence, we obtain that, for any e > 0, 



E f T \\Z r \\ 2 p{X r ) 4 / q+£ dr<oo. 
Jo 



In the next proposition we find an adapted process smaller than Y. This 
process will give us a lower bound on the explosion rate of Y on the blow-up 
set {£ = oo}. 

Proposition 5 (Lower bound on Y). We define the following process: 

1 \ 



Then for all t>0, St<Yt 



j(t-t At) + l/& 



Proof. Denote by at the quantity 

1 \ l / q 



a t 



,q(r-TAt) + 

if t < t and at = £ on {t > r}. The process a solves the equation 

a t = a T -J^ al +q l[ 0tT ](r)dr. 

Note that H$ = ~E{at\J-t)- Thanks to Jensen's inequality, 

S t <E^^H T - f t B}+n [0>T] (r)dr 

and the comparison theorem (Corollary 4.4.2 in [6]) achieves the proof. □ 

We now want to find a lower bound for p(Xt/\ T ) 2 ^ q St when t goes to oo 
on {£ = oo}. 

Lemma 3. Let p(x) denote the distance of x G D to the boundary dD 
and t be the exit time from D of the diffusion X. Let the conditions (L), 
(B) and (E) hold. Then there exist two positive constants C\ and Ci which 
depend on D, q, a and b such that for all x £ D, 

IN 1/91 



d < p(x) 2 / q E x 



<C 2 
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Proof. Recall that if x G D, ~P x (t > 0) = 1 and 



If r < h, then sup tg r ^ \X t — x\ > p(x). Therefore, we can apply Theorem 
4.2.1, page 87 of [23] to obtain 

P x (t <h)<P x ( sup \X t -x\> p(x)) < Kl e K2h e- K2p W 2 / h . 
\te[o,h) / 

We apply this inequality with h = 1 jy q and y > 1 : 



E, 



^) < i + K ie K ^y q e- K ^ 2 y q dy 

r+oo 

<l + K ie K2 / e- K ^ yq dy 



qp(x) 2 /l 

Since — 1 + 1/q > — 1, we deduce 



f>K2 r+oo 

qp(xy/ q J p (x)2 



For the other inequality remark that 
and P x (t <l/y q )> P x {X 1 / yq D). We just have to find a lower bound to 

, r l/p(x) 2 /l r+oo _ A„. 

p{xf/ q / P x {X 1/yq iD)dy= P x {X p{x?uq i D)—. 

Jo J 1 u 

Let A x be the set (R d \D) n B(x,2p(x)) which is not empty, and VolfAj.) 
denotes the volume of the set A x . Using the Aronson estimates of [22], we 
have, for u > 1, 

1 \<*/2 

>^3 ' 



2vruV(^) 2 
x e W (-K 4 u q p(x) 2 - K ^J p X {iy ) dy 
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Jt 



d/2 



A, \ 2iTuip{x)' 
/ I \d/2 /■ 



*M-K^ P ( x y) p{x)d 

because u > 1. Thus, 

2/9 ^ V*{X x ,i,tD)dy 

Vol(A a 



exp(-K4U 9 /9(x) ) 



u 2+dq/2 



p{x) 



with 



/I \ d / 2 

* = . 

The integral has a lower bound because the open set D is bounded and the 
dominated convergence theorem shows that 

f + °° t , ^2^ du f + °° du 

We have supposed that dD € C 2 (which was important in the proof of The- 
orem 6). Therefore, the curvature is continuous on dD which is compact; so 
the curvature is bounded. There exists r > such that each point y € dD 
lies on the boundary of a ball with radius r and this ball is contained in the 
complementary of D (see Figure 1). 

Instead of calculating the volume of (M. d \D) H B(x, 2p(x)), our problem 
is reduced to the following: we find the volume of the intersection of two 
balls in that case; see Figure 2. 

If r < x < 3r, the volume is equal to 

pa(x) rPix) 

V(x) = Cr d {sm6) d d9 + C2 d {x-r) d {sm9) d d6, 
Jo Jo 

where C is the volume of the unit ball in R d_1 , 



a(x) = Arccos 



x 2 + r 2 — 4(x — r) 
2xr 
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boundary of 
ball of radius r 




Fig. 1. 



and 



(5{x) = Arccos 



/ hx — 3r 



Now we must prove that 



(x—r) d 



> K\. We split V(x) in two parts. For the 

first part, C2 d (x — r) d ^^ x \s\nO) d d6 ', the result is clear because if r < x < 
2r, 

< Arccos(7/8) < < -. 

3 

For the second part, Cr d Jq' (sin 9) d d6, if r < x < ^^- r, 
x 2 + r 2 — 4(x — r) 2 



2xr 



€[0,1] 



a(x) € [0,vr/2] 



and we use the fact that sin is an increasing function on [0,7r/2[; so 

p*tf m < «w (i - (" 2 + r2 2 ;f "° 2 V ^ 
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„2 L „2 a!^ „\2\d/2 



<2 d / 2 a(x)(l 
3 \ d / 2 



2xr 



( 3 \ 1 

— a(x)(x - r) . 
\xr J 



Therefore, if r < x < 4+v ^ r, 



F(ar) 
(a; — r)' 



>C2 d 



Arccos(7/8) 



(sin0) d d0>O. 



This proves that > K and therefore, p(x) 2 E x (±) > d = K 5 K. This 

achieves the proof in the uniformly elliptic case. □ 

Remark 5. If the diffusion matrix is degenerate, the result on the lower 
bound may be false. Suppose that a = and b is bounded by k. If the exit 
time r is smaller than l/y q , 



k_ . r^y" 
~y q 



sup \X t — x\> p(x) 
[o,i/y q ] 



ri/y ft 
> / \b(X r )\dr> sup / b(X r )dr 
Jo [o,i/y] Jo 

and thus, 

and the limit, as p{x) goes to zero, is zero. 

From the inequality (12), we already know that there exists a constant C 
such that 



Vt > Y f < 



C 



(piXtAr)) 2 /* 



Now we prove Proposition 1: 



Proof of Proposition 1. From Proposition 5, we work on the process 
3 = (St)t> . For alH>0, 

1/9 



p y q (X tAT )E t = p^(X tAT )E^ 
= p 2 / q (X tAT )E^ 
+ p 2 ^(X tAr )E^ 



q(r - t A t) + 

t q M<oo 



l + q£«(T-T At) 
1 



[(r-rAt) 



4=oo 
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The first term in the right-hand side is nonnegative. Let t = t — t At: it is 
the first exit time of the diffusion starting at Xt/\ T . Hence, 

r/ i 

p 2 /*(X tAT )E^ ' 



[(r-rA t) 



1\V9 



E^V /a (^tAr)E 



T 



>Ci - 



1\ 1/9 



E^(l{=oo), 



where C\ is the lower bound of Lemma 3. Thus, we obtain 



p 2 /«(X Mr )~ t > d 



and we deduce the announced result. □ 



E' 



1-5=00 ) 



3. Continuity. Recall that we have constructed a couple of processes 
(Y, Z) which satisfy for all rj > and all < t < T, Y t > and 



YtAT v — YTATr, 



TAtv, 



1At„ 



(y r ) 1+9 dr 



TATn 



Zj> dBy. . 



tAT„ 



Moreover, on the set {t > r}, Y t = ^, Z t = and liminft_> +00 Yi Ar > £ a.s. 
We now want to prove the converse inequality, namely, limsup t _ ++00 Y t A T < £ 
a.s. Remark that we just have to show this estimate on the set {£ < +00}. 

3.1. Existence of the limit. We first prove that the limit of Yt^ T , as t 
goes to +00, exists a.s. In the proof we will distinguish the two cases: £ is 
greater than a positive constant and £ is nonnegative. 

3.1.1. The case where £ is bounded away from zero. We can show that 
(XtA T )t>o has a limit as t — > +00 by using Ito's formula applied to the process 
l/(Y n ) q . We prove the following result: 

Proposition 6. Let the conditions (B) and (E) hold. Suppose there 
exists a real a > such that £ > a > 0, P-a.s. Then 

-l/q 



(20) 



YtAT 



( ^ 



0<t, 



where is a nonnegative supermartingale such that on the set {t > t}, 
*t = 0. 



Proof. From Proposition 5, for every n £ N* and every < t, 



y t n > ~™ = E' 



L\g(r-rAt) + (l/(^ An))« 
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Since £ > a, we have 
Z r t l > E** 



> a 



g(r-rAt) + (l/a)' 
1 N 1 /? 



l/q- 



>aE :Ft 



1 



1 + qrofi 



l/q 



1 + ga^E-^ (r) 



Therefore, 
(21) 



Vt > < ^ < 1(1 + gaflS*(T)) < +oo 



because the conditions (B) and (E) hold, which implies, in particular, that 
r G L^fi). Thus, for all 4 > 0, {Y t n )- q belongs to L 1 ^). We want to apply 
the Ito formula to the semi-martingale Y n with the function < x i— > x _<? . 
But we just have that a.s. for all t > 0, Y/ 1 > 0. For e > 0, we define a C 2 
function f £ : R — > R such that on R + , 

'■<*> = (iTi)' 

Note that for a fixed x € R+, (f e ( x ))e>o is increasing and the limit is equal 
to f(x) = x~ q . By the Ito formula, for all < t < T, 



(22) 



fs^Xtkr) = feO^TAr 

TAt 
tAr 



TAi 



tAr 



f' £ (Y r n )Z r n dB r 



TAt 



tAr 



f' £ (Y r n )(Y r n ) l+q dr 



fl AT 

-\\ f' £ '(Y?)\\Z?\\ 2 dr 

JtAr 

i-TAt 

E*/ 6 p<? AT ) -E* / f' £ (Y r n )(Y r n ) 1+q dr 

JtAr 



TAt 



tAr 



f' £ '(Y r n )\\Z?fdr. 



1 \ 2+q 



x + e 



Now for x > 0, /i(x)x 1+< ? = -g(^) 1+9 , thus, -q < f' £ (x)x l+q < 0, and 

£'(*) = g(l + . 
Thereby, a.s. and in 2^(0) for all < t < T, 

rTAr 

limE^ / f^Y r n )(Y r n ) 1+q dr = -qE^(T Ar-t At) 

JtAr 

From (21), we have that a.s. and in L 1 (f2) 

limE^/ e (y T V)=E^— I 



(Y^ArY 
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and 



lim/ e (lft T ) 



£-♦0 



(Y t n A r) q 



For the last term in (22), we use the monotone convergence theorem and 
hence, we have proved that, for all < t < T, 

«b^'mb +?E * (TAT - 4AT) 



1(1 + 1) E r, l T " r M 



TAr \\yn\\2 



Let T go to +oo: 



1 



E 



1 



tAr 



+ qW t {r -tAr) 



(yn)2+. 



dr. 



(23) 



n||2 



Mr (W+* 



dr. 



Let n>m. Since (Aii>fAm, we obtain, for all < i, 
1 1 

< 



1 1 



E 



(fAm)9 ((An)« 
E" 



q(q + l) f^ f r \\ZT\\ 2 



ds - E Tt 



\Z 



n\\2 



tAr 



{Y n )q+ 2 



ds). 



Now 



+ 1) 



E 



m\\2 



tAr (n m ) 9+2 



ds - W?* 



\\Z. 



n\\2 



< 



E 



(£Am)i ((An)! 



V 



tAr (Y?)^ 
1 



1 



For a fixed t > 0, the sequences (E^ ^/^^ )n>i and Qyi • )q )n>i converge a.s. 

and in L 1 (dominated convergence theorem). Then the sequence (E^/^ r ^y^J +li ds) n >i 
converges a.s. and in L 1 and we denote by $ the limit 

|| Z n„2 



= lim 



On the set {£ > r}, <I>t = a.s. and we have 



{Y n )q+ 2 



ds. 



g(g+i) E ^ r _E 



n||2 



tAr 



ds 
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< qE^ir -fAr)+E ft [ - 1 . 

\(£ An) 



Thus, 



For r <t, 



* t <gE^(r) + — . 



E^ / '' ' " rfs > E^E^ / '' \ " ds, 
$ r >E^'-$t. 

We deduce that ($t)o<t is a nonnegative supermartingale. Now for all n € N* , 

Fix i > 0. Taking the limit as n — > +oo, we deduce 

9E- Ft (r-tAr)+E- Ft f^- > ) - $ t . 



(*W) 9 

This achieves the proof of Proposition 6. □ 

$ being a nonnegative supermartingale, the limit of $tAr as t goes to +oo 
exists P-a.s. and this limit <£ r - is finite P-a.s. The L 1 -bounded martingale 
E^'^) converges a.s. to as t goes to +oo, then the limit of Y t ^ T as 
t — > +oo exists and is equal to 

lim Y tAr - 



t^+oo IAT (l/^-$ T -)V?' 

If we were able to prove that <3? is continuous (or $ T - is zero a.s.), we would 
have shown that Y is a continuous process. 



3.1.2. The case £ nonnegative. Now we just assume that £ > 0. We can- 
not apply the same arguments because Y n may to equal to zero with positive 
probability, which implies, in particular, that (Y t n )~ q L 1 (f2). We will ap- 
proach Y n in the following way. We define for n > 1 and m > 1, £ n > m by 

f' m = KAn)V-. 
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This random variable is in I? and is greater or equal to 1/m a.s. The BSDE 
(3) with £ n > m as terminal condition has a unique solution (y"'™,^'" 1 ). It 
is immediate that if m < m! and n<n' , then 

As for the sequence Y n , we can define Y m as the limit when n grows to +00 
of Y n ' m . That limit Y m is greater than Y = lim ?w+oc Y n . But for m<m! for 
0<t <T, 



^v m - ^v m = « - « - / [iXr n n q+l - or™ )* 1 ] *■ 

/•TAt _ _ , 

- / [Z r n ' m - Z r n ' m ] d£ r 

J!At 

< y™^ _ y»V _ [ TAT [Z?>™ - Z?> m '\ dB r 

JthT 

and taking the conditional expectation given 

Letting first T — > +00 and then m' — > +00 in the last estimate leads to 

n < y n ' m _ Y n < — 

Therefore, P-a.s., 

su P |y^ r -y t Arl<-. 

t>o m 
Since for each m > 0, (y^ T ) f >o has a limit on the left at +00, so does Y". 

3.2. Continuity ofY. We know now that 
(24) liminfy Ar = lim Y tAT > £ 

and on the set {t > r}, y = £. In this part we give sufficient conditions to 
ensure that the process Y is continuous, that is, 

lim Y Mr = £. 

t— >+oo 

It suffices to prove the result on the set {£ < 00}. In the rest of this section, 
we will suppose that (Al) and (A2) hold, and P(£ < 00) > => ^ <9D. 
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3.2.1. A first step. In the first section we have proved the following es- 
timate: 

(12) P-a.s. Vt>0 \Y t \< C 

(/WatJ) ,q 

where p is the distance to the boundary of D. The constant C depends on 
q, D and the bound on b and cr in (B). Here we want to construct another 
estimate which depends also on the function g. Our result is the following: 

Proposition 7. Suppose that the boundary of D belongs to C 3 . If U 
is an open set such that U n = and U fl dD ^ 0, then there exists a 
constant C = C(U,g,q,b,o~,D) and an open set Djj such that D C Djj and 
if Pu denotes the distance to the boundary of Djj, we have 

(25) P-a.s. Vn G N, Vt > Y? < ° 

{pu{XtAr)r /q 

Recall that r is always the first exit time from D. 

Proof. We suppose that the set = {g = +00} is not equal to dD. 
Hence, if we define for all e > the set 

F £ = {yedD; dist(y, F^) < e}, 

there exists e' > such that F e i 7^ dD. Moreover, if U is an open subset of 
R d such that F^^U = and U ndD ^ 0, there exists < e < e' such that 
UndDcdD\F e . 

Recall that D is a bounded open set of M. d with a boundary dD £ C 3 . 
Thus, there exists r > such that on T r = {y G R d ; dist(y, <9D) < r}, the 
signed distance d 

( dist(x,dD) = p(x), if x G D, 
\-dist(x,dD), if x£R d \D, 

belongs to C 3 (r r ). Moreover, for all y G r. r , there exists a unique x G dD 
such that y = x — d(y)Tt(x), where Jt(x) is the outward normal vector at 
the point x G dD. We have \\y — x\\ = \d(y) \ = dist(y, dD). The result can be 
found in [9]. 

We take a function i/j e :M d —> [0,1] such that ip £ is of class C 2 (R d ) and 
ip £ = on F e and ^ E = 1 on 3D \ i^e • With this function we define the set 
D e as follows: 

D £ = D U {y G E d ; 3x G dD, 3u G [0,r/2[ s.t. y = » + uip £ (x)lt(x)}, 

where rf (x) always denotes the outward normal vector at the point x G dD. 
We can easily prove that D £ is included in Dur r , and that 

F £ C dD £ and dD \ F 2e C D £ . 



d(x) 
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If dD £ C 3 , then the boundary of D £ is of class C 2 , and from our construc- 
tion, the distance to the boundary of D £ , denoted by p £ , is also a C 2 function 
on the set T r . Moreover, if y G dD \ F2 £ , then p £ {y) = r/2 > 0. 

Now the proof of (25) is similar to the proof of Theorem 6. Let <3? be 
a C°°(]R d ) function such that $ = 1 on D \ T r and $ = on T r/2 . For all 

< r] < r/2 and C > 0, we define a function * = $> v € C 2 (M d ;]R + ) such that 
on D, 

T T C C 



" [(l-^ft + i^ + T?] 2 /? [0„]2/9- 

The constant i? e is the supremum of /0 £ on D £ and we can easily see that 
Re < sup p + r/2 and that 6^ > p £ for all 77 > 0. Remark also that ^ is of class 
C 2 on D. We apply the ltd formula to ^{Xt/\ T ) and by the same arguments 
as in the proof of Theorem 6, we can choose the constant C [depending only 
on D, on q and on the bound of b and a in (B)], such that for all < t < T, 

i-TAt 

(-XW) = ^(Xtat) - / VV(X r )a(X r ) dB r 

JtAr 

pTAt rTAi 



P l AT rt AT 

/ m{X r ) 1+q dr + / C/ r dr; 

JtAr JtAr 



HAt JtAr 

with [/ a nonnegative adapted process. The constant C must satisfy (15), 
that is, 

20 1 / 2 \ 

C q + —lCVe„)b - - - + 1 I IkWJI 2 + -2 Trace (ctct* Z> 2 0„) > 0. 



Moreover, on {t > t}, 



because on dD, <3? = 0, and on F £ , p £ = 0; 

because on <9L>, < p £ < r/2. Recall that for all n € N, (Y n ,Z n ) is the 
solution of the BSDE (3) with terminal time r and terminal data g An. On 
the compact set dD \F £ , by (A2), the function g is bounded by a constant 
-K" = K £ . We choose C > and < r/ < r/2 such that 

c c 

> n and , —. > K. 



ifli ~ (ri + r/2) 2 li 

We can take C > Kr 2 l q satisfying (15), and rj < r/2 A C q / 2 /n q / 2 . Note that 
C does not depend on ij. Therefore, if we define, for t > 0, 

Y t = V(X tAT ) and Z t = V*(Jft)<7(X t )l t<TJ 
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the process (Y, Z) satisfies P-a.s., for all < t < T, 

„ _ pTAr _ pTAr pTAr _ 

Y tAT = Ytat ~ / Y r 1+ i dr + U r dr- Z r dB r , 

JtAr JtAr JtAr 

U being a nonnegative process, and on the set {t > r}: Yt > g{X T ) An. From 
the comparison theorem (Corollary 4.4.2 in [6]) for solutions of a BSDE, we 
obtain 

P-a.s. Vt >0F f "< * v (X tAr ) < ° 

{p e (x tAT )y/i 

Since this inequality holds for all n, we have proved the proposition. □ 

The main interest of Proposition 7 is that if h £ Cq(U) (h has a compact 
support included in U) with UnF^ = and UndD ^ 0, then the sequence 
h(X. AT )Y n is bounded in L°°([0, +oo[xf2): there exists a constant K = Kjj 
such that 

P-a.s. Vt > h(X tAr )Y t n < K. 
We can also prove the following: 

Proposition 8. For all u > 1, there exists a constant K = Kum > 
such that 



E 



r\\z t f P 4 j q+ \x t )dt<K. 

Jo 



Proof. Using Proposition 7, the proof is the same as the proof of Propo- 
sition 4. □ 

3.2.2. Continuity: the conclusion. Recall that Foo = {g = +00} n dD is 
a closed set, that U is an bounded open set such that U H F^ = and 
UndD^0. 

Now we take a function tp : M. d — > M + of class C 2 and with a compact 
support included in U. For /3 > 0, we apply the Ito formula to the process 



(26) 



E(e-^( 5 An)(X r MX r )) 



-0E [ T e-P r <p{X r )Y?dr + E [ T e- pr v{X r )Y?\Y?\ q dr 

JtAr JtAr 

+ E T e-P' r Y r n £ip(X r )dr + E f e~ Pr Z?Vy{X r )a(X r ) , 

JtAr JtAr 
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where C is defined by (7). In (26), every term, except maybe the last one, is 
well defined because f3 > 0, ip is a C 2 function with compact support, and 
Y n is bounded by n. Now using the Cauchy-Schwarz inequality, we obtain, 
for all r] > 1, 



E 



e- pr \Z?- V<p(X r )a(X r )\dr 



(27) 



< 



E 



i T \\zn 2 pT\x r )dr 

Jo 



1/2 



E 



2/3r -4/g-»? 

Pu 



(X r )\\S7<p(X r )(r(X r 



• dr 



1/2 



We already know that a is bounded on D. The support of Vtp is included 
in U. In our previous construction of Djj, we have U D D G D\j and on 
U H D, pu > r/2 > 0. Therefore, Pj/^ ^Vc^ is a continuous and bounded 
function. With Proposition 8, we deduce 



(28) 



E 



[ T e~P r \Z? ■ V<p(X r )a(X r )\dr < K. 
Jo 



It is important to remark that the constant K does not depend on n. 

We want to pass to the limit when n — > +oo in (26). With the monotone 
convergence theorem, we obtain, for all < t, 

lim E( e - 8T ( J An)(X>(Ir))=E(e^(I>(I T )); 

n— >+oo 



lim E 

n— >+oo 



e- pr cp(X r )Y r n dr = E / e'^ ip{X r )Y r dr; 



tAr 



-/3r, 



lim E 

n— >+oo 



tAr 



n\l+q 



dr = K 



tAr 



tAr 



-? r <p(X r )(Y r ) 1+ idr. 



The support of the function C<p is included in U. Therefore, from Proposi- 
tion 7, Y n dp{X) is a.s. bounded. From the dominated convergence theorem, 
we deduce that, for all t > 0, 

lim E f T e~ l3r Y r n £ip(X r )dr = E f e- pr Y r &p(X r ) dr. 

ra->+oc J tAr J tAr 

The last term in (26) is equal to 



E 



Mr 



e -/3r Z n. W{X r )G{X r )dr 



E 



tAr 



e -Pr p 2/ g +v/2 {Xr)z n . p - 2 /^/ 2 ( Xr )V^>(X r ) dr. 



From Proposition 8, the sequence p 2 jj q+r ^ 2 (X)Z n l T> . is bounded in L 2 ([0, 
+oo[xfi) for all r/ > 1. Therefore, after extraction of a suitable subsequence, 
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which we omit as an abuse of notation, p^/ 9+r ^ 2 (X)Z n l T> . converges weakly 
in L 2 ([0, +oo[xfi). The process 



p. -2/q- V /2 
Pu 



T>- 



is in L 2 ([0, +oo[xfi) because a is bounded, and from our construction, 
Py 2 ^ 9 (X)V <p(X)l T> . is also bounded. Using Proposition 3, we obtain 



lim E / e- pr Z?-V<p(X r )<j(X r )dr = E eT^ Z r ■ Vcp(X r )a(X r ) dr. 

n->+oo J tAr J tAr 

Finally, letting n — > +00 in (26), we have, for all < t, 
E(e-P T g(X T )<p(X T )) 

= E(e-^Y tA MX tAT )) 

(29) 

-PE f e-P r <p(X r )Y r dr + E [ e~P r y{X r ){Y r ) 1+q dr 

JtAr JtAr 

+ E f e-P r Y r &p{X r )dr + E [ e'^ Z r ■ V<p(X r )a{X r ) dr. 

JtAr JtAr 

From Proposition 7, we know that 

E ( e~ Pr tp(X r )Y r dr <K, 
Jo 

E f e-P r ^{X r ){Y r ) 1+q dr + E f ' e -^Y r \Op{X r )\dr < K. 



For the last term, using the Cauchy-Schwarz inequality [see (27)] and Propo- 
sition 4, we obtain 

E f T e~ Pr \Z r ■ V<p(X r )a(X r )\dr < K. 



Therefore, when t goes to +00 in the equation (29), we obtain, using 
Fatou's lemma, 



(30) 



E(e-^g(X T )ip(X T ))= lim E(e~^ tAT ~>Y tAT ip(XtAr)) 

t—f+OC 



> E 



e~%(X r )( lim Y tAr 



But recall that we already know that 



(24) 



lim Y tAT >g{X T )- 

t— >+oo 



Hence, the inequality in (30) is in fact an equality, that is, 



E(e- l3r g(X T MX T ))=E 



e-f T ip(X T )[ lim Y tAT 
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And using again (24), we conclude that 

lim Y tAr = g(X T ), P-a.s. on {g{X T ) < oo}. 

t — >+oo 

4. Minimal solution. In the third section we have constructed a process 
(Y, Z) which satisfies the conditions (Dl) and (D2) of Definition 2. We will 
prove now that, if this process is a solution of the BSDE (3), that is, if it 
satisfies also the condition (D3), then it is the minimal nonnegative solution. 



Theorem 7. Let the conditions (L), (B) and (E) hold and let (Y,Z) be 
a nonnegative solution of the BSDE (3) (solution in the sense of Definition 
2). Then P-a.s. for all t > 0, 

Y t > Y t . 



Proof. Recall that is the first exit time of D \ r„ and we have, for 
all < t < T, 



TAt 



<At„ 



(Y r ) 1+q dr 



TAt 



(At. 



Z r dB r 



For ?jGN*, (Y n ,Z n ) is the solution (in the sense of Definition 1) of the 
BSDE (3) with terminal data (An, We compare Y with Y n : 



YtAr v - ^tAT, — ^TAr^ - ^TAt,, 
fTAr,, 



TAt„ 



tAr„ 



\Y r ) l+q -(Y r n ) 1+q dr 



(At. 



(Z r -Z™.)dB r 



Ytat v - Ytat v 

cTAt„ 



TAt,-, 



tAr„ 



a?(Y r -Y r n )dr 



(At„ 



(Z r — Z™) dB r 



where the process is defined by 



Q 



(V — (Y n \ 1+ i 



Y —Y n 

1 r 1 r 



a? = (l + q)(Y?)*, 



if Yr=Y". 



a n is a nonnegative process and we have a linear BSDE whose solution is 

rTAr„ 



(31) Y tATv -Y t n AT =K- 



(Xta, ,-l? Ar Jexp 



ctldr 



tAr„ 
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From the hypothesis of the theorem, Y is nonnegative and Y n is bounded 
by n. Indeed, on the set {t > r}, Y™ = £ A n < n and for all < t < T, 

i-TAt i-TAt 

Y t l T = F T " AT - / (Y r n )^ dr - / Z? dB r 

JtAr JtAr 



TAt 
tAr 



<Y£ Ar - I Z?dB r 



thus, 



y™ T < Y? - [ T Z? dB r = £ A n - [ T Z? dB r <n- f Z r r l dB r . 

JtAr JtAr JtAr 

Taking the conditional expectation, we deduce that, for all t > 0, Y t n < n. 

We now pass to the limit in (31) first as r/ — > 0, then as T — > +oo and with 
the Fatou lemma, we obtain for all t > 0, 

Therefore, Y is greater than Y n for all n G N* and thus greater than Y. □ 
Moreover, we obtain the following result: 

Proposition 9. There exists a constant C (the same constant as in 
Theorem 6) such that P-a.s., for all t>0, 

C 



Y t < 



tAr) 



Proof. For all sufficiently small rj > 0, we denote by p v the distance 
from the boundary of D v = D\T ri , that is, 

D v = {xeD,p(x)>rj}. 
If x € Drj, p(x) — rj < prj(x) < p(x). We consider the first exit time 

T v = in{{t>0,X t <£D v }. 
From Theorem 6 we deduce 

Vt>0 FtAr,<-27^ nVKY \ W 

Pr,' q (X tATv ) P < q {XtAT v ) -V 

The constant C which appears in the previous inequality may depend on rj. 
In the proof of Theorem 6 we use the fact that there exists p > such that 
on r„, the signed distance function is of class C 2 . But if rj < fi, it is also 
true that p v is of class C 2 on r„. So in the proof of the theorem we can use 
the same function (p and the same bound R for and p. Moreover, on T^, 
\Vp\ = 1 and D 2 p depends only on the curvature of dD. Therefore, we can 
choose a constant C independent of r] if r] < p. 

To conclude, let rj — > and we obtain the desired inequality. □ 
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5. Viscosity solution of the associated elliptic PDE. Recall that D is a 
bounded open subset of M. d with a C 3 boundary. For all x £ D, {X x ;t > 0} 
is the solution of the SDE (4): 

(4) Xf = x+ f b{X x )dr + f a(X x )dB r for t > 0. 

Jo Jo 

The functions b and cr are defined on M. d , with values respectively in R d 
and M dxrf , and such that b and cr are continuous on M rf and satisfy the 
conditions (M), (L) and (B). For each x G D, we define the stopping time 
t x = mf{t > 0, Xf £ D}. We assume that 

(32) P(t x < oo) = 1 for all xeD, 
that the set of singular points 

(CI) r = {x G 3D; P(t x > 0) > 0} is empty, 

and that for some (3 > and all x G D, 

(C2) Ee^ < oo. 

Let us recall the following result (cf. Proposition 5.2. in [18]): 

Proposition 10. Under the conditions (CI) and (C2), the mapping 
x i — ► t x is a.s. continuous on D. 

Let g : dD — > M + be a continuous function and for all n G N, we define 
9n = g An. Hence, g n is a continuous function. For all n G N, from Remark 4, 
{(Y t x,n , Z^ ,n );t > 0} is the unique solution (in the sense of Definition 1) of 
the BSDE (3) 

(33) Y t x ' n = g n (XD- r Y r *' n \Y r x > n \Ur- H Z? n dB r . 

JtAr x Jt/\T X 

We denote by u n the function defined on D by 

u n {x) = Y x ' n . 

For h G C(dD,M), we consider the elliptic PDE (6) with boundary con- 
dition h: 

— Cv + v\v\ q = on D; 

v = h on dD. 

The following definition can be found in [1] and [2] (or [5] and [18] for v 
continuous). If v is a function defined on D, we denote by v* (resp. v*) the 
upper- (resp. lower-) semicontinuous envelope of v. for all x G D, 

v*{x) = limsup v{x ) and v*{x) = liminf _v(x'). 
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Definition 3 (Viscosity solution). 

• v : D — > R is called a viscosity subsolution of (6) if < +oo on D and if 
for all <fi £ C 2 (R d ), whenever x £ D is a point of local maximum of v* — <f>, 

-£(f)(x)+v*{x)\v*{x)\ q <0 ifx£L>; 

min(-£</)(x) + t>*(x)|v*(x)|V*(x) - h(x)) <0 if x € 3D. 

• v : D — > R is called a viscosity supersolution of (6) if t>* > — oo on D and if 
for all (f) £ C 2 (R c( ), whenever x G D is a point of local minimum of v* — (j>, 

-C(j)(x) +v*(x)\v*(x)\ q > ifxGD; 

max(— C(j)(x) + v*(x)\v*(x)\ q ,v(x) — h(x)) > if x £ 

• u : Z) — > R is called a viscosity solution of (6) if it is both a viscosity sub- 
and supersolution. 

Let us recall the following result (cf. Theorem 5.3. in [18]): 

Theorem 8. Under the assumptions (M), (L), (B), (CI) and (C2) ; 
since g An is continuous on dD, u n is continuous on D and it is a viscosity 
solution of the elliptic PDE (6) with boundary data g An. 

Remark 6. Since g A n is continuous on dD, from Theorem 3.3 in [5], 
it follows that u n is the unique continuous viscosity solution of the PDE (6) 
with terminal data g An. 

From now on we add the uniformly elliptic condition: there exists a con- 
stant a > such that, for all x £ R d , 

(E) aa*(x) >ald. 

With this assumption, (CI) and (C2) hold if (B) is true. In the previous 
sections we have constructed a process {(Y t x , Zf);t > 0} which is a solution 
of the BSDE (3) with terminal data g(X* x ) (in the sense of Definition 2). 
Y x is the limit of Y x > n : for all t > 0, 

(16) Y?= lim Y? n . 

n— >+oo 

If we define 

u(x)±Y*, 

then u is the limit of the sequence u n . Thus, u is nonnegative. Since u is the 
supremum of continuous functions u n , u is lower-semicontinuous on D and 
satisfies 

C 

(34) Vx £ D u(x) < 
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Recall that p is the distance from the boundary dD and C is a constant which 
does not depend on g. Moreover, u(x) = g(x) on dD. Since g is not bounded 
on dD, we cannot apply Theorem 8. Moreover, the condition v* < +00 in 
Definition 3 cannot be satisfied on D. Therefore, we change the definition 
of a solution. 

Definition 4 (Unbounded viscosity solution). We say that v is a vis- 
cosity solution of the PDE 

— Cv + v \v \ 9 = on D, 

(6) 

v = g on dD, 

with unbounded terminal data g if v is a viscosity solution on D in the sense 
of Definition 3 and if 

g(x) < lim u*(V) < lim v*(x')<g(x). 

x'— >x x'— >x 
x'£D,x£dD x'eD,x£dD 

Remark that this definition implies that v* < +00 and > —00 on D. 
5.1. u is a viscosity solution. 

Lemma 4. The function u is a viscosity solution of the PDE (6) on D. 

Proof. We will use the half-relaxed upper- and lower-limit of the se- 
quence of functions u n : 

u(x) = limsupit n (x') and u(x) = liminf u n (x'). 



X'— *x 



x —*x 



Since {u n } is a nondecreasing sequence of continuous functions, we have 
\/x £ D u(x) = u*(x) = u(x) < u*(x) = u(x). 

We fix r] > and we prove that on D = D\{x £ D, p(x) < 77} , u is a viscosity 
solution. We already know that 

~ C 
Vx € D u(x)<—— r . 

Recall that continuous viscosity solution and from the Lemma 6.1 of 

[5], we deduce that u is a viscosity solution of (6) on D and this holds for 
all i] > 0. Therefore, the lemma is proved. □ 



Since u n is a nondecreasing sequence of C°(D) functions, we have, 
(35) \/x£dD liminf u(x ) > g(x) = u(x). 

x 1 — *x, x'&D 

Hence, u* is a supersolution of (6) because u* > g on dD. 
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Lemma 5. The solution u satisfies the boundary condition, that is, 

lim u*{x) < g(x) = u(x). 



X — >x 
x'eD, i63D 



Proof. We already know that 

liminf u(x') > g(x) = u{x). 



X'— >x 

x'eD.xedD 



So we just have to prove the converse inequality on the set {g < +00}. If U 
is an open set such that U D -Foo = and U D dD 7^ 0, there exists an open 
set Djj and a constant C\j such that, for all n € N, 

Cu 



P-a.s. 



vt > y^ ,n < 



Recall that /?[/ is the distance to the boundary of Djj. From the proof of 
Proposition 7, the choice of the set D\j and of the constant C\j does not 
depend on x S D. 

We write again equation (26): 

E(e-^( 5 An)(X*>(X*J) 

= u n (x)cp(x) 

(36) 

- /?E r e-^ifiX^Y^ dr + E F" e~ pr -<p(X?)Y?> n \Y?> n \ q dr 
Jo Jo 

+ E r e-^Y^&piX?) dr + E f" e~ pr Z* >" • V^(X r x )(j( X r x ) dr. 
Jo Jo 

The function ip : M d -> M + is of class C 2 and has a compact support included 
in U. The constant /? is positive. From Proposition 7, there exists a constant 
K[/ such that, for all n £ N, 



E 



<i^t/E 



Moreover, using (27), we have 



E 



e- pr \Z^ n V^{X^)a{X^)\ dr 



< 



E 



71,71 II 2 4/g+77 



p^"(X?)dr 



1/2 



E 



- 2 ^p- 4/9 -"(^)||V^(^>(^)|| 2 dr 



1/2 
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Recall that p^^ q is a continuous and bounded function on D. With 

Proposition 8, we obtain that 



E 



[ TX e- pr \Z*> n ■ Vip(X*)a(X*)\dr < KuE f' e~ Pr dr. 
Jo Jo 



Since x^ t x is a continuous function on D and since t x = if x E dD, 
we have 

lim E f Tx ' e _/3r dr = 0. 



X — >X 







If x G and if (x m ) me pj is a sequence of elements of D which converges to 
x, we replace in (36) n by m and x by x m and we take the limit asm-* +oo. 
We obtain, by Fatou's lemma, 

limsupu m (x m )v9(x m ) = limsupE(e- /3T ^ {g A m)(X^ )ip(X*™J) 

m— >+oo m— >+oo 

< e( limsupfe-^- (g A m)(X*™ M**™ )f 

\ m— »+oo 

because gy? is a bounded function. By continuity of x *— > and of gip, we 
have 

limsupn m (x m )93(x) = limsupn m (x m )(^(x m ) < g(x)(p(x). 

m-^+oo m— >+oo 

Finally, on {g < oo}, we have 

limsup u*(x')<g(x). 

x f — *x 

x'eD, xedD 

With inequality (35), this achieves the proof of the lemma. □ 

5.2. Some regularity results on u. We want to prove now that u is con- 
tinuous on D. Here it seems to be necessary to assume the condition (E). 

First we prove that, under stronger assumptions on b and a, u belongs to 
C°(D;R + )nC 2 (D;R+). 

Proposition 11. Recall that b and a satisfy always (L), (B), (E) and 
dD £ C 3 . We assume that b and a belong to C 1 (D). Then u is in C 2 (-D;R+). 

In order to prove this result, we need the following lemma: 

Lemma 6. The assumptions of Proposition 11 hold. We consider a con- 
tinuous function h : dD — > IR. Suppose that (Y, Z) is the solution (in the sense 
of Definition 1) of the BSDE 

Y t = h(X T )- [ T Y r \Y r \ q dr- f Z r dB r . 
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Then there exists a function v : D — > R + of class C°(D) n C 2 (D) such that 

Vt>0 Y t = v(X tAT ) and Z t = Vv(X tAT )a(X tAT )l t<T . 
Moreover, v is solution of the PDE (6) with boundary condition h. 

Proof. Since h is continuous and since b and o belongs to C 1 (D), from 
the Theorem 15.18 in [9], there exists a unique solution v 6 C°(D) (1 C 2 (D) 
of the PDE (6) (see also [15]). We prove that, for all t > 0, Y t = v{X tAr ) 
and Zt = V v{X tAT )a{X tAT )lt<T ■ We want to apply the Ito formula to the 
process v(X). But we just have v € C 2 (D) and we do not know if we can 
define a function v E C 2 (M rf ) such that v = v on D. 

We will use some arguments of the proof of Theorem 15.18 in [9] . We define 
a sequence {h m } of functions such that {h m } approximates h uniformly on 
dD and h m 6 C 2 ' 7 (IT). From Theor em 15.10 in [9], there exists a function 
v m such that u m solves the Dirichlet problem (6) with condition h m on 
the boundary and v m € C 2n (D). We apply the Ito formula to the process 
v m (X): for all 0<t <T, 

rTAr 

v m (X tAT ) = v m {X TAT ) - \ (Cv m )(X r ) dr 

JtAr 

TAt 

Vv m (X r )a(X r ) dB r 

rl At 



= v m (X T AT)- v m {Xr)\v m {X r )\ q dr 

JtAT 

[TAt 

- I Vvm{X r )a(X T )l r<r dB r . 

JtAT 

We denote by (Y m ,Z m ) the solution of the BSDE (3) with terminal data 
h m (X T ) £ L°°(f2). Uniqueness of solution of this BSDE implies 

Vi>0 Y t m = v m (X tAT ) and Z™ = Vv m (X t )a(X t )l t<T . 

From (C2) and Remark 4, we know that there exists a constant C such that, 
for all m € N, 



0<Yj n = v m {x) <E 



sup e^|F t m | 2 
*e[o,r] 



<(M[e pT \h m {X T ) 



Sinc6 h m converges uniformly to h on dD^ hm 

is a bounded sequence in 
L°°{dD). Therefore, the sequence {v m } is uniformly bounded on D. 

From Theorems 6.1, 13.1 and 15.3 in [9], the sequence {v m } converges 
uniformly on compact subsets of D, together with its first and second deriva- 
tives, to the function v. 



42 A. POPIER 

Since {v m } is uniformly bounded on D and converges to v, for all < t < 

T, 

(38) lim Y t m = lim v m (X tAr ) = v(X tAr ), 

m— >+oo m— >+oo 
/•TAr /-TAt 

lim / v m (X r )\v m (X r )\Ur = / W (X r )|w(X r )|«dr. 

nw+oo J tAT J tAr 

Using Ito's formula and the Burkholder-Davis-Gundy inequality, we obtain, 
for a constant c independent of m, 



E 



sup \Y t m -Y t \ z + I \\ZT - Z t \\ z dt 

0<t<T 



<cE\h m (X T )-h(X T )\ 2 . 



Therefore, with (38), we conclude that a.s. Yt = v{Xt/\ T ) for all t > 0. More- 
over, there exists a constant K such that 

E [ T \\Vv m {X r )a{X r )l r<T \\ 2 dr = E [ \\Z™\\ 2 dr < K < +oo; 
Jo Jo 

and with (37) and (38), for all < t < T, 

i-TAt rTAr 

lim / Vv m {X r )a{X r )l r<T dB r = I Z r dB r . 

m^+oo J tAr J tAr 

Let KZ be a compact subset of D. Since the first derivatives of v m converge 
uniformly on KZ, from the dominated convergence theorem, we deduce 

rTAr 

lim / \\(Vv m (X r )a(X r )-Vv(X r )a{X r ))l r<T l,c(X r )\\ 2 dr = 0. 

Therefore, for all compact subset KZ of D, P-a.s., 

Z t l K {X t ) = Vv(X t )a(X t )l t<T l K (X t ). 

If {KZ m } is an increasing sequence of compact subsets of D such that \J m KZ m = 
D, for all m, 



E 



f T \\Vv{X t )(j{X t )l t<T l Km {X t )\\ 2 dt<E [ \\Z t \\ 2 dt < +oo 
Jo Jo 



and since r > t implies Xt G D, with the monotone convergence theorem, we 
deduce 



E 

Then 

E 



f \\Vv{X t )a{X t )l t<T \\ 2 dt<oo. 
Jo 



r \\Vv(X t )a{X t )l t<T - Z t \\ 2 dt = 0. 
Jo 



This achieves the proof of the proposition. □ 
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From Lemma 6, we can deduce that u n belongs to C°(D) n C 2 (D) if a 
and b belong to C 1 (D). 

Proof of Proposition 1 1 . We fix rj > and we consider the set D v = 
D\{x G D, p{x) < r]} for all n > 0. From Lemma 6, u n G C 2 (D) and satisfies 
—Cu n + «^ +9 = on D. And on Z^, it n is bounded by C /rj 2 / q . Therefore, 
from Theorem 15.5 in [9], we obtain that S7u n is bounded on Z^n- The 
sequence is bounded in C 1 (D2 V ), thus the limit u is continuous on D2 V , that 
is, u is continuous on D. Moreover, we already know that u is continuous 
on the boundary. Therefore, we deduce that u belongs to C°(D,]R+). 

Now if we consider the PDE, 



from Theorem 15.18 in [9], the equation has a regular solution v G C°(D r? ) n 
C 2 (D V ). But this solution is also a continuous viscosity solution. Since u is 
now a continuous viscosity solution of the same PDE, from the comparison 
result in [5], we deduce that v = u, that is, u G C 2 (L>^). Hence, u belongs to 



Now we want to prove that u is continuous on D without the regularity 
conditions on b and a of the Proposition 11. We just assume that (M), (L), 
(E) and (B) hold. 

Proposition 12. The viscosity solution u is continuous on D and is 
locally Holder continuous on D. 

Proof. We will show that for all open sets D' C D such that D' C D, 
there exists < a < 1 such that the sequence of functions u n is bounded in 
the space C a (D'). C a (D') is the set of functions v such that 



Since u n converges to u, we deduce that u belongs to C a {D') and thus is 
continuous on D. 

In order to prove that u n is a bounded sequence in C a (D'), we will con- 
struct a sequence v m which will belong to C a {D') and such that there exists 
a constant K such that, for all m G N, ||u m || Q < K. Let b m and a m be two 
sequences of functions such that: 

1. b rn and a m belong to C 1 (D) and b m and a m are bounded in L°°(D); 

2. b m (resp. a m ) converges to b (resp. a), uniformly on D; 



Cv — v\v\ q = 



on D, 



v = u 



on dD 



C 2 (D). □ 
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3. a m satisfies the condition (E). 

Let v m be the unique solution in C^(-D) H C 2 (D) (see Lemma 6 or [9]) of 
the equation 

- C m v m + v m \v m \ q = on D, 

(6) 

v m = g An on <9D, 

where £ m is the operator: 

Vx£R d C m ip{x) = \Tra,ce(<j m a* m {x)D 2 v{x)) + b m (x)V<p(x). 
For xeD, let X x > m be the solution of the SDE 

Vt > Xf ' m = x + /* 6 m (X*' m ) dr + f a m (X^ m ) dB r , 
Jo Jo 

r m is the first exit time from D of the diffusion X x ' m , (Y x > n < m , Z x > n > m ) is the 
solution of the BSDE: 

Y t x ' n ' m = (g An)(X x ^ n )+ J Y x ' n ' m \Y x > n > m \ q dr - J Z x ^ m dB r . 

From classical results on the SDE, X x,m converges to X x solution of the 
SDE (4) and the process (Y x > n > m , Z x ' n ' m ) converges to (Y x ' n ,Z x ' n ) solution 
of the BSDE (33) (see Proposition 4.4 in [6]). 
From Lemma 6, we have 

v m {x) =Y x ' n ' m and u n {x) = Y^ ,n . 

Therefore, v m converges to u n . Moreover, we know that v m is a bounded 
sequence in L°°(D). 

Let D' be a open subset of D such that D' C D. We apply Theorem 8.24 
in [9]. The function v m is the solution of 

£v m = v m \vm\ q = g G L°°. 
Therefore, there exists a real < a < 1 and a constant K such that 

IbmlU < -K"||Vm||L°°- 

The constants depend on the ellipticity constant of a m , on the bound on b m 
and a m in L°° and on the distance between D' and dD. We deduce that u n 
belongs to C a (D') and the norm ||n n || a is bounded w.r.t. n E N. 
Finally, u belongs to C a (D'). □ 
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5.3. Minimal viscosity solution. We prove the following: 

Theorem 9. If v is another nonnegative viscosity solution of the PDE 
(6) {in the sense of Definition 4), and if v* > g on dD, then u<v on D. 

Proof. We show that for all n G N*, u n <v*. The proof is the same as 
the proof of Theorem 3.3 in [5]. We fix n, we assume that there exists z G D 
such that 5 = u n (z) — v*(z) > and we will find a contradiction. The main 
tool is Theorem 3.2 in [5]. □ 

6. Other generators /. We have considered the generator f{y) = —y\y\ q . 
The main properties of this function are it is nonincreasing and allows the 
explosion at time r [see (5) for the definition of this stopping time]. But we 
can also consider more general generators. Let be a nonincreasing 

function of class C 1 , such that there exists q > 0, k > s.t. 

(9) Vy > f(y) < -Ky 1+ «. 

The BSDE (10) has a unique solution if £ satisfies the condition (H6). From 
Remark 4, if £ G L°°, then (H6) holds. We also assume that /(0) = 0; thus, 
if £ > 0, then Y t > for all t > 0. 

First of all, the conclusion of Theorem 6 holds: there exists a constant C 
such that for every solution (Y,Z) of the BSDE (10), 

V420 |y ' |s w£))^ 

Indeed, with the notation of the proof of Theorem 6, equality (13) becomes 

TAt 



pi AT 

*(X tAr ) = *(At At ) + / f(*(X r )) dr 

JtAr 



TAt 

W(X r )a(X r )dB r 

tAr 
TAt 

[V^(X r )b(X r ) 

ItAT 

+ 1 Tvace{aa*(X r )D 2 ^(X r )) + f(V(X r ))] dr. 
With assumption (9), equation (14) becomes 

-(Vtf)ft- ^Trace(o-o-*Z9 2 *) - /(*) 



10 

KC q + —(X79)b 

q 

+ ||crV0|| 2 + - Trace(crcr*D 5 

q\q J q 
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and we can choose C such that the right-hand side is nonnegative. The rest 
of the proof remains the same. 

Now we suppose that £ is a nonnegative, JP" T -measurable random variable 
such that P(£ = +00) > 0. As in Section 2, we construct a process (Y,Z) 
satisfying the conditions (Dl) and (D2) of Definition 2: (Y, Z) is the limit of 
the sequence of solutions (Y n , Z n ) of the BSDE (10) with terminal condition 
£ A n. 

For the continuity of Y [condition (D3) of Definition 2], Section 3.2 re- 
mains unchanged, if we have already proved that the limit of Yt AT when t 
goes to +00 exists. 

We define the following function F on K_: 



With (9) and since / is of class C 1 and nonincreasing, F is a positive, de- 
creasing and convex function such that lim^o F(y) = +00, and lim^+oo F(y) = 
0. Moreover, for e > 0, F £ is defined by F £ (y) = F(y + e). 

Now if a is a constant such that £ > a > 0, for all n G N, all t > 0, Y t n > 
a.s. We can apply the Ito formula to F £ (Y n ): for all < t < T, 



Now F" > and < F' £ (y)f{y) = < 1. Let T go to +00: 

F £ (Y t n Ar ) < E^'F £ (£ Art) + E^*(t - t Ar) < F £ {a) +E^*(t - t A r) 



Recall that from (C2), r G L x (0). Hence, for all t > 0, sup n F(y t n ) belongs 
to L 1 . With (39) and the same ideas as in the proof of Proposition 6, we 
deduce that, for t > 0, 



where $ is a nonnegative supermartingale. 

Finally, if / is a nonincreasing and C 1 function with /(0) = 0, such that 
(9) holds, and if £ is a nonnegative, ^--measurable random variable such 
that P(£ = +00) > 0, then the BSDE (10) has a minimal solution (in the 
sense of Definition 2). And the associated PDE (11) has also a minimal 
viscosity solution. 

The only thing which we cannot describe just with inequality (9), is the 
behavior of Y on the set {£ = +00} (see Proposition 1). 




(39) 





< F(a) +E >t (r-t Ar). 



F{Y tAT ) = E * (F(0 + r - t A r) - $ 
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